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ABSTRACT 

We investigate stellar microlensing of the collimated gamma-ray burst 
afterglows. A spherical afterglow appears on the sky as a superluminally 
expanding thin ring ("ring-like" image), which is maximally amplified as it 
crosses the lens. We find that the image of the collimated afterglow becomes 
quite uniform ("disk-like" image) after the jet break time (after the Lorentz 
factor of the jet drops below the inverse of the jet opening angle). Consequently, 
the amplification peak in the light curve after the break time is lower and 
broader. Therefore detailed monitoring of the amplification history will be able 
to test whether the afterglows are jets or not, i.e., "disk-like" or not, if the 
lensing occurs after the break time. We also show that some proper motion 
and polarization is expected, peaking around the maximum amplification. The 
simultaneous detection of the proper motion and the polarization will strengthen 
that the brightening of the light curve is due to microlensing. 



Subject headings: gamma rays: bursts — gravitational lensing 



1. INTRODUCTION 



Whether the gamma-ray bursts (GRBs) are collimated or not is one of the most 
important questions in GRBs, since the jet configuration is crucial for almost all aspects 
of GRBs, such as the energetics, statistics, and central engine models of GRBs. There are 
several suggestions for the collimation of GRBs. When the Lorentz factor of the ejecta 
drops below the inverse of the opening angle, the ejecta begins to spread sideways, and we 
expect a break in the light curve of the afterglow (Rhoads 1999). Such breaks have been 
observed in several GRBs, such as GRB 990510 (Stanek et al. 1999; Harrison et al. 1999) 
and GRB 991216 (Halpern et al. 2000). The rapid decline rate (e.g., GRB 980519) or the 
unusually slow decline rate of GRB 980425 also suggests the collimation of GRBs (Halpern 
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et al. 1999; Sari, Piran & Halpern 1999; Nakamura 1999). However, in order to establish 
that the afterglows are jets, other observations are indispensable, such as early afterglow 
observations at radio frequencies (Frail et al. 2000). Microlensing of the afterglows has 
a potential to be one of such observations. The emission region of the afterglow seen by 
an external observer occupies a size of ~ 10 17 cm on the sky which is comparable to the 
lensing zone of a solar mass lens located at the cosmological distance (Loeb & Perna 1998). 
Therefore microlensing can resolve the afterglow, and hence we can expect that some 
features of jets can be obtained by the observations of the microlensed afterglows. 

On the other hands, afterglows will be useful to constrain the cosmological density 
parameter of stellar mass objects through microlensing. As a result of the relativistic 
motion, the emission region appears to expand superluminally, so that the microlensing 
time is only ~ 1 day (Loeb & Perna 1998). This is much shorter than that of common 
sources, ^ 10 years (Gould 1995). The microlensing probability of a source at a redshift 
z s ~ 3 is ~ 0.3fi* (Press & Gunn 1973; Gould 1995; Koopmans & Wambsganss 2000). 
Since the known luminous stars amount to ^ 4 x 10~ 3 (Fukugita, Hogan & Peebles 
1998), roughly 1 in ~ 10 3 GRBs would be within the lensing zone. If some fraction of the 
dark matter is made of massive compact halo objects (MACHOs), as Galactic microlensing 
searches suggest (Alcock et al. 2000), we can expect higher probability. Remarkably, the 
microlensed afterglow may have been already observed in the light curve of GRB 000301C 
as the achromatic bump (Garnavich, Loeb & Stanek 2000; Gaudi, Granot & Loeb 2001; 
but see also Panaitescu 2001). Since new GRB satellites such as HETE-2 and SWIFT wi\\ 
provide many opportunities to monitor afterglows frequently, it is important to study the 
jet effects on the microlensing of the GRB afterglows. 

In §2 and 3, we will calculate the equal arrival time surfaces, and the surface brightness 
distribution of the afterglow image, respectively, taking into account of the jet effects. In §4, 
we will calculate the microlensed light curve. In §5, we will investigate the proper motion 
and the polarization induced by microlensing. §6 is devoted to summary and discussions. 



Let E, p, 6i and c s be the burst energy, the density of the ambient gas, the initial 
opening angle, and the sound speed of the jet, respectively. The Lorentz factor of the shock 
front approximately evolves as 



2. EQUAL ARRIVAL TIME SURFACE 




r < r b , 
r > r b , 



(1) 
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where r b = (75E /Mpcl) 1 ^ = 3.9 x lO 18 ^ 3 ^ 1/3 (v^c s /cT 2 / 3 cm, r r = (8/75y/ 3 r b and 
T b = 2c s /5c9i (Rhoads 1999). Here E 52 is the burst energy in units of 10 52 ergs, and ri\ 
is the ambient gas density in cm -3 . We assume that the material is uniformly distributed 
across the jet at any r, the emission comes just behind the shock front, and r 3> 1. 
Such approximations are sufficient and convenient to understand the key features before 
performing more detail calculations (Rhoads 1999; Panaitescu & Meszaros 1999; Moderski, 
Sikora & Bulik 2000; Huang, Dai & Lu 2000), although the future detailed analysis may 
demand the realistic model (Granot, Piran & Sari 1999; Granot & Loeb 2001; Gaudi, 
Granot & Loeb 2001). 

The arrival time T of a photon emitted at an angle 9 from the line of sight is given by 
T = t — rfi/c where t — J dr/c\/l — T~ 2 and /i = cos 9. This relation with equation ([1]) can 
be written as 
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r > r b , 
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where T b = r b /8cT 2 b = 0.35EH\ 1/3 (v / 3c s /c)- 8 / 3 (^/0.01) 2 days. Equation <g) gives 
the equal arrival time surface at T, and some examples are shown in Figure 1. In 
the power law regime, T ^ T b , the apparent size r^ max of the jet is determined 
by the equal arrival time surface (Sari 1998; Panaitescu & Meszaros 1998). The 
distance of the equal arrival time surface from the line of sight is given by 
r± = r\J\ — /i 2 ~ r^2(l — /i) with equation (H), which takes the maximum value 

of r ±tmax = (T/5T b ) 5 / 8 r b /T b = 1.2 x 10 17 E 5 1 2 /8 ^r 1/8 (^/0.01)- 1 / 4 (T/lday) 5 / 8 cm at 
r / r b = (T ' /hT b ) l l A . On the other hand, in the exponential regime, T ^ T b , the apparent size 
of the jet is determined by the opening angle of the jet, 9 b = 9i + c s t co /r, where t co = J dr/cT 
is the comoving time in the jet frame. With equation (|l]), we find 
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r < r b , 
r > r b . 



(3) 



If the direction of the jet axis is (9, 4>) = (#o, 0) with 9 = being the line of sight, the edge 
of the jet is given by 9 2 + 9q — 299 q cos0 = 9 b . In Figure 1, the trajectories of jet edges for 
9 = and 9 = 9. t on the = plane are plotted. We can see that there is an intermediate 
regime where the evolution of the Lorentz factor is power law while the edge of the jet 
determines the apparent size (Rhoads 1999; Panaitescu & Meszaros 1999). 
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3. 



SURFACE BRIGHTNESS DISTRIBUTION 



Let us consider a system that has an isotropic luminosity dL u in its rest frame and is 
moving with Lorentz factor 7 in a direction (9, 0) with a solid angle dfl = dfidcj). Emitted 
photon with frequency v is blueshifted to frequency — (3fi), where (3 = y/1 — 7~ 2 . 

Using the Lorentz transformation, the observed flux at distance D and frequency v is given 
by dF v = dL vl (i_p^ j 'An D 2r y 3 (1 — /3fi) 3 . A jet, with a luminosity L u in its local frame, is a 
collection of such systems, so that we have dL u = L u dQ/ir9 2 with equation (^). Now we 
assume that the luminosity in the local frame is L v oc T a v b . Then, the flux per unit area, 
i.e., the surface brightness S(r±, 0; T) := d 2 F u /r±dr±d(fi can be obtained by 



where we use dQ/r^dr^d<f) = 2dfi/dr 2 L and the fact that the Lorentz factor of the shock 
front T is higher than that of the material behind it 7 by a factor of \/2 (Blandford & Mckee 
1976). With equations fll]), (@) and @, we can find that equation (Q) is a function of r, 
noting r_i_ ~ rJ2(l — p) and 1 — /3/i ~ 1 — (3 + 1 — /i~ 1/T 2 + 1 — p. Since we can obtain r 



for given r± by r± ~ 2(1 — /x) with equation @, the surface brightness can be calculated 
as a function of r±. Note that there are two solutions of r for each r±, one from the front of 
the equal arrival time surface and the other from its back. As time goes, the emission from 
the back does not contribute due to the jet geometry, as we can see from Figure 1. 

The assumption L v oc Y a v h is valid for frequencies far from the typical synchrotron 
frequency v m (Sari, Piran & Narayan 1998; Sari, Piran & Halpern 1999). For simplicity, 
we neglect scattering, self-absorption, and electron cooling. Then, the luminosity is 
proportional to the total number of swept-up electrons N e oc T~ 2 (Rhoads 1999) times 
the radiation power from each electron P oc 7 2 I? 2 oc T 4 . Since the typical frequency is 
u m oc -B7 2 oc T 3 , the luminosity at the typical frequency is L Vm oc At v <^ u m , we 
have L u = L i/m (z//z/ m ) 1 / 3 oc T -2 // 1 / 3 . At v ^> v m , if electrons are accelerated to a power law 
distribution with index p, we have L v = L Um {u / Um)'^ 1 ^ 2 oc r( 3 P~ 5 )/ 2 z/~( p-1 )/ 2 . Therefore, 
we find a = —2, b = 1/3 at v <^ u m , and a = (3p — 5)/2, b = —{p — l)/2 at v v m . We 
shall adopt p = 2.5. 

In Figure 2, we show the surface brightness distribution as a function of the distance 
from the center for 9 = (the jet axis coincides with the line of sight). In the power 
law regime T/Tb <^ 1, the surface brightness is brighter near the edge and dimmer at the 
center, i.e., "ring-like" (Waxman 1997; Sari 1998; Panaitescu & Meszaros 1998). However, 
in the exponential regime T/Tb ^ 1, the surface brightness becomes nearly constant, i.e., 
"disk-like". The qualitative feature does not depend on 9q. In the intermediate regime, a 
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part of the ring will be missing due to the jet edge when 6*0 ^ 0. 

4. MICROLENSED LIGHT CURVE 

We now consider a point mass lens of mass M that is located at L from the 
source center on the sky (see Mao & Loeb 2001 for binary lenses). The Einstein 
radius on the source plane is given by R E = [(4:GM/c 2 )(D s Di s /Di)} 1 ^ 2 , with D[, D s 
and Dis being the angular diameter distance to the lens, to the source, and from 
the lens to the source, respectively (e.g., Schneider, Ehlers & Falco 1992). A point 
lens is described by three parameters, the Einstein radius, r E := R E / (r b /T b ) = 

0.45 (M/M ) 1/2 [( J D s A s /A)/10 28 cm] 1/2 E 52 1/3 r2i /3 (VZcs/cf'* (^/O.Ol) -1 , the impact 
parameter, I := L/(r b /T b ), and the azimuthal angle of the lens position, <pi, with respect to 
the source center on the sky. Here we measure the distance from the line of sight in units 
of r b /T b . 

The observed flux can be calculated numerically by F v = J A(r±, <f>)S(r±, 0; T)rj_dr^d<fi, 
where A(r±, 4>) = (uf + 2)/uiyJuf + 4 and u x = [L 2 + r\ - 2Lr ± cos(</> - </>/)] 1/2 / R E . The 
unlensed flux can be obtained by putting A = 1. In Figure 3, the lensed flux of an afterglow, 
the unlensed flux, and the lensed flux while retaining the initial surface brightness (i.e., 
retaining the "ring-like" image) are plotted by solid, dashed, and dotted lines, respectively. 
Only the exponential regime T ^ T b is depicted, where the slope of the unlensed light 
curve is — p = —2.5 at v ^> u m , and —1/3 at v v m (Sari, Piran & Halpern 1999). At 
early times, the temporal profiles of the lensed and unlensed flux have the same shape but 
different amplitudes, since the source can be regarded as pointlike. The offset between the 
lensed and unlensed flux is a function of l/r E - The maximum amplification occurs around 
when the edge of the image crosses the lens, r_|_ m(M . ~ L = lr b /Y b . Here we set I = 1 so that 
the the maximum amplification occurs at T ~ 10T&, with the image being "disk-like". 

In the upper panel of Figure 3, we consider the case v 3> v m . At the maximum 
amplification, an amplification peak would appear if the image were "ring-like" , as in the 
power law regime[| (Loeb & Perna 1998). However, since the image is "disk-like", the 
amplification peak becomes lower and broader, and apparently disappears for v ^> u m . 
While for v u m , as shown in the lower panel of Figure 3, the low amplification peak does 
exit in the light curve even if the image is "disk-like" . This is because the decline rate of 
the unlensed flux at v ^> v m is too steep (— p = —2.5) for the time derivative of the flux to 



1 This statement may depend on the more detail profile of the ring (Granot et al. 1999). 
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be positive, but the slope at v <C v m is very mild (—1/3). Note that the amplification factor 
is nearly the same in both cases. In Figure 4, the amplification factor (while retaining the 
initial surface brightness) is plotted by solid (dotted) lines. The amplification factor for the 
uniform ring with a fractional width 10%, which reproduces the bump in the light curve of 
GRB 000301C (Garnavich, Loeb & Stanek 2000), is also plotted by dashed lines. Errorbars 
of ±0.1 magnitude are also shown for reference. As analyzed by Gaudi & Loeb (2001) and 
Gaudi, Granot & Loeb (2001), a future monitoring campaign of a lensed afterglow will 
be able to reconstruct the radial structure of the afterglow image, although the quality of 
the observational data for GRB 000301C (about ~ ±0.1 magnitude) may not be sufficient 
for this purpose. If the reconstructed image is "disk-like" after the break time, it will be 
strengthened that the afterglows are jets. 

In the intermediate regime, where the evolution of the Lorentz factor is power law but 
the edge of the jet determines the image size, the amplification depends on the lens position 
<pi when 6q ^ 0. If the lens is located on the side where the ring is cut off by the jet edge, 
the amplification peak is small, while, if the lens is on the other side, the peak is relatively 
large. 



5. PROPER MOTION AND POLARIZATION 

Even when the multiple images formed by microlensing can not be resolved, 
the centroid of the combined image is expected to move (Hosokawa et al. 1993, 
H0g, Novikov & Polnarev 1995, Walker 1995, Miyamoto & Yoshii 1995, Paczyhski 
1998, Mao & Witt 1998). The order of the maximum displacement is estimated by 
R E /D S ~ 10- n (M/M Q )^ 2 ~ 2(M/M ) 1 / 2 /ias with R E ~ lO 17 (M/M ) 1 / 2 cm and 
D s ~ 10 28 cm. This proper motion might be measured by the upcoming missions, such 
as the Space Interferometry Mission (SIM; see [http:/ /sim.jpl.nasa.govp, with positional 



accuracy down to ~ 1/xas (Paczyhski 1998). To neglect the proper motion due to the jet 
edge effect (Sari 1999), we set here 9o = 0. We can calculate the light centroid numerically 
by weighting the image positions with brightness (Walker 1995, Mao & Witt 1998). The 
upper panel of Figure 5 shows the proper motion as a function of time at v ^> v m . Initially, 
since the source is pointlike, the initial displacement is a function of 1/te (see Figure 1 in 
Walker 1995). As the source expands, the centroid moves toward the lens due to the finite 
source effect (Mao & Witt 1998). The centroid goes through the origin in coincidence with 
the maximum amplification. The upper panel of Figure 6 shows the maximum magnitude 
of the displacement as a function of the normalized impact parameter 1/te for the uniform 
surface brightness (thick lines), the initial surface brightness (dotted lines) and the uniform 
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ring with a fractional width 10% (dashed lines). Since the surface brightness is retained, the 
displacement depends only on 1/te- The displacement is smaller for more uniform surface 
brightness but is always larger than the initial displacement. 

We also expect some polarization in the microlensed afterglows (Loeb & Perna 
1998). Here, following Sari (1999), we shall consider the favorite conditions in which 
the polarization at each point in the image is n = (p + l)/(p+ 7/3) ~ 72% toward 
the source center. To neglect the polarization due to the jet edge effect (Ghisellini & 
Lazzati 1999; Sari 1999), we set here 6q — 0. The net polarization can be calculated by 
averaging Hocos(20) with A(r±, <j))S(r±, <p; T). The lower panel of Figure 5 shows the net 
polarization as a function of time at v ^> v m . The polarization has a maximum toward the 
lens in coincidence with the maximum amplification. The lower panel of Figure 6 shows 
the maximum polarization as a function of the normalized impact parameter 1/te for the 
uniform surface brightness (thick lines), the initial surface brightness (dotted lines) and the 
uniform ring with a fractional width 10% (dashed lines). Even in the exponential regime, 
the polarization due to microlensing is comparable to the polarization which is expected 
around the jet break time when 6q ^ (see Figure 4 in Sari 1999). 

6. SUMMARY AND DISCUSSIONS 

We have investigated stellar microlensing of collimated GRB afterglows taking into 
account of the jet effects. Using the analytical expressions for the evolution of the jet, we 
have first calculated the surface brightness distribution of the jet. We find that the image 
is "disk-like" after the jet break time, rather than "ring-like" before the break time. We 
have further analyzed the microlensing signal in the light curve of the afterglow observed 
at frequencies far below and far above the typical synchrotron frequency u m , where the 
radio is below v m for about 1 month,Q and the optical is above v m after about 1 day. If the 
edge of the image crosses the lens before the break time, the microlensing signal appears as 
an achromatic amplification peak in the light curve (Loeb & Perna 1998). The peak seen 
in GRB 000301C is believed to be a microlensing event, occurred before the break time, 
in the "ring-like" regime (Garnavich, Loeb & Stanek 2000; Gaudi, Granot & Loeb 2001; 
but see also Panaitescu 2001). We find that, after the break time, the amplification peak 
becomes lower and broader because of the "disk-like" image. Since detailed monitoring of 
the amplification history will be able to reconstruct the afterglow image (Gaudi & Loeb 
2001; Gaudi, Granot & Loeb 2001), it could be tested whether the afterglows are jets or 



2 The radio scintillation is suppressed at v ;> 10 GHz (Goodman 1997). 
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not, i.e., the reconstructed image is "disk-like" or not after the break time. We should take 
care that the amplification peak apparently disappears for v ^> v m . At the break time, the 
peak amplitude depends on the lens position. 

Microlensing also induces the proper motion and the polarization of the afterglow in 
coincidence with the maximum amplification. The magnitude of the proper motion is order 
of the Einstein radius of the lens, which might be detected by upcoming missions, such as 
SIM, with positional accuracy down to ~ 1/xas. From the proper motion we can estimate 
the apparent size of the afterglow explicitly. The maximum polarization due to microlensing 
is comparable to the polarization which is expected around the jet break time (Sari 1999), 
even if the image is "disk-like" . The simultaneous detection of the proper motion and the 
polarization will strengthen that the brightening of the light curve is due to microlensing, 
although an initial constant positive offset in the amplification history (Loeb & Perna 1998) 
and a specific level of chromaticity to the amplification history (Granot & Loeb 2001) are 
unique features of the lensed afterglows. 

We would like to thank the anonymous referee for useful comments to improve our 
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Fig. 1. — Equal arrival time surfaces, i.e., the surfaces from which photons from a thin shell 
arrives at the same time T to the observer far on the right side, are shown by solid lines. We 
adopt 6i = 0.01 and c s = c/a/3. The upper and lower panels are the same, but the scale of 
the vertical axis is different. The arrival times are T/T5 = 10 -2 , 10 _1 , 1, 10, 10 2 as the surface 
gets larger. The Lorentz factor of the emitting shell obeys equation ([!]), in which the shell 
decelerates exponentially at r > r b . The trajectories of the jet edges for 9 = and 9 = 61 
are also plotted with dotted and dashed lines, respectively, with 8 = being the line of 
sight. No photon comes from the outside of these trajectories. 
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Fig. 2. — The surface brightness as a function of the distance from the center at the times 
T/Tf, = 10 -2 , 10 _1 , 1, 10, 10 2 is shown for frequencies above and below the typical synchrotron 
frequency with solid and dotted lines, respectively. Each surface brightness is normalized 
by the average surface brightness. The relative brightness at the center increases with the 
time T. The image is "ring-like" at T £ T b and "disk-like" at T £ T b . The divergence of the 
brightness at r± = r± jmax is an artifact of the assumption that the radiation comes from a 
two-dimensional shell. 
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Fig. 3. — The lensed flux of a GRB afterglow, the unlensed flux, and the lensed flux while 
retaining the initial surface brightness at all the time as a function of the time T are plotted 
by solid, dashed and dotted lines, respectively. The origin of the vertical axis has no meaning. 
We set te — 1, I — 1, Oq — and c s = c/ VS. The upper panel shows the flux at frequencies 
above the typical synchrotron frequency, and the lower panel is at frequencies below the 
typical synchrotron frequency. In the upper panel, the boxed region is expanded. 
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Fig. 4. — The amplification factor (while retaining the initial surface brightness) as a function 
of the time T is plotted by solid (dotted) lines. The amplification factor for the uniform ring 
with a fractional width 10% is also plotted by dashed lines. Errorbars of ±0.1 magnitude 
are also shown for reference. We set te — 1, I — 1, 0q — and c s = c/^/~3. The upper panel 
shows the flux at frequencies above the typical synchrotron frequency, and the lower panel 
is at frequencies below the typical synchrotron frequency. 
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Fig. 5. — In the upper panel, the light centroid position of a GRB afterglow on the sky in 
units of the ratio of the Einstein radius Re to the source distance D s as a function of the 
time T is shown for three cases (r E ,l) = (0.1,0.03), (1,0.2) and (1, 1). The horizontal axis 
has the same meaning as in the lower panel. In the lower panel, the polarization of a GRB 
afterglow as a function of the observed time T is shown for (r E , I) = (0.1, 0.03), (1, 0.2) and 
(1,1). The polarization scales linearly with the polarization at each point of the image Il , 
and we assume here n = (p + l)/(p + 7/3) ~ 72%. In both panels, frequencies are above 
the typical synchrotron frequency. 
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Fig. 6. — In the upper panel, the maximum magnitude of the light centroid displacement 
as a function of the normalized impact parameter 1/te is shown for the uniform surface 
brightness (thick lines), the initial surface brightness (dotted lines) and the uniform ring 
with a fractional width 10% (dashed lines). The horizontal axis has the same meaning as 
in the lower panel. In the lower panel, the maximum polarization as a function of the 
normalized impact parameter l/r E is shown for the uniform surface brightness (thick lines), 
the initial surface brightness (dotted lines) and the uniform ring with a fractional width 10% 
(dashed lines). We assume here n = (p + l)/(p+ 7/3) ~ 72%. In both panel, the surface 
brightness is retained, so that the displacement and the polarization depend only on 1/te- 



